Mon. Not. R. Astron. Soc. 000, HUTU (20091 Printed 23 February 2009 (MN M£X style file v2.2) 



Formation of the Radio Profile Components of the Crab 
Pulsar 



o 
o 

(N 

X> 
<D 
IX, 

m 

(N 
P? 

43 

9^ 
6 

CO 

03 



> 
ON 

00 

m 

(N 
O 
On 
O 



S. A. Petrova * 

Institute of Radio Astronomy, NAS of Ukraine, 4> Chervonopraporna Str., 61002 Kharkov, Ukraine 



Received. 



ABSTRACT 

The induced Compton scattering of radio emission off the particles of the ultrarela- 
tivistic electron-positron plasma in the open field line tube of a pulsar is considered. 
We examine the scattering of a bright narrow radio beam into the background over a 
wide solid angle and specifically study the scattering in the transverse regime, which 
holds in a moderately strong magnetic field and gives rise to the scattered component 
nearly antiparallel to the streaming velocity of the scattering particles. Making use of 
the angular distribution of the scattered intensity and taking into account the effect of 
rotational aberration in the scattering region, we simulate the profiles of the backscat- 
tered components as applied to the Crab pulsar. It is suggested that the interpulse 
(IP), the high-frequency interpulse (IP') and the pair of the so-called high-frequency 
components (HFC1 and HFC2) result from the backward scattering of the main pulse 
(MP), precursor (PR) and the low-frequency component (LFC), respectively. The 
components of the high-frequency profiles, the IP' and HFCs, are interpreted for the 
first time. The HFC1 and HFC2 are argued to be a single component split by the 
rotational aberration close to the light cylinder. It is demonstrated that the observed 
spectral and polarization properties of the profile components of the Crab pulsar as 
well as the giant pulse phenomenon outside of the MP can be explained in terms of 
our model. 
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1 INTRODUCTION 

The radio emission pattern of the Crab pulsar is unique 
in its complexity. There are seven emission components lo- 
cated at different longitudinal regions throughout the whole 
pulse period and characterized by distinct sp ectra and polar- 
ization (e.g. iMoffett fe Hankinsf[l996l . 119991 ). At the lowest 
frequencies, the radio profile of the Crab pulsar consists of 
three components: the main pulse (MP) , the precursor (PR) 
~ 15° ahead, and the interpulse (IP) ~ 150° behind the MP. 
Such a structure is also characteristi c of the profiles of some 
other pulsars (e.g., PSR B 1055-52. iMcCulloch et all 119761 
and PSR B1822-09, iFowler. Wright fc Morris! Il98ll ). How- 
ever in the Crab pulsar the components have exception- 
ally steep spectra, with the power law indices qpr = —5, 
qip =4.1 and «mp = —3, and vanish at the frequencies 
~ 0.6, 3 and 8 GHz, respectively. In the range 1-3 GHz, the 
Crab profile contains the so-called low-frequency component 
(LFC) ~ 36° in advance of the MP. At still higher frequen- 
cies, there appear the high-frequency interpulse (IP') ~ 10° 
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earlier in pulse longitude than its low-frequency analogue 
(the IP) and a pair of the so-called high-frequency compo- 
nents (the HFC1 and HFC2), which lag the IP' by ~ 70° 
and 130°, respectively. The IP', HFC1 and HFC2 have flat 
spectra and dominate the pulse profile at the frequencies 
> 8 GHz. 

The radio polarization studies of the Crab pul- 
sar ([C ampbell, Hcilcs & Rankin 1970; Manchester 197ll: 
Manchester. Huguenin fc Taylor! Il972l; IMoffett fc Hankinsl 



19991 : iKarastergiou. Jessner fc Wielebinskil 12004! ) have re- 
vealed moderate linear polarization of the MP, IP and LFC 
(~ 15%, 25% and 40%, respectively), high polarization of 
the HFC 1,2 and complete linear polarization of the PR and 
IP'. In the MP, PR and IP, the position angle of linear po- 
larization is nearly constant across the components and is 
nearly the same for all of them. In the rest of the com- 
ponents and between them, the position angle noticeably 
changes with pulse longitude. In the LFC and IP', the char- 
acteristic values of the position angle differ from that of the 
MP by ~ 30° and 90°, respectively, whereas in the HFC1 
and HFC2 the position angle is shifted by ~ 90° from that 
of the LFC. A straightforward interpretation of the posi- 
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tion angle swing in terms of the primord ial rotating vec- 
tor model dRadhakrishnan fc Cooke 19691) faces di fficulties 
ijMoffett fc Hankindll999l : lKarastergiou et al.ll2004l ). 

In spite of substantial distinctions in the spectra and 
polarization, most of the components exhibit similar fluctu- 
ation properties. The Crab pulsar is known as a source of 
giant pulses - the occasional strong radio pulses with inten- 
sities ranging from a few dozen to a few thousand times the 
average (e.g.. iLundgren et al.|[l995h . The gia nt pulse activ- 
ity is characteristic of the MP, IP and IP ' (|Sallmen et alj 
1 19991 ; ICordes et a"ill2004l ; iKong et alJl20Qgft. and it has also 
been found in the HFC1 and HFC2 (fjessner et aill2005l ). In 
these components, the giant pulses may be present simulta- 
neousl y, but the rate of the i r occurrence is sub s tantially dif- 
ferent dSallmen et al.lll999l ; ICordes et ai1l2004 IPopov et all 
l200d : lKong et al.ll2008l) . Furthermore, the temporal and fre- 
quency structure of the giant emission in t he MP and IP' is 
essentially distinct (|Hankins fc Eilek|[2007h . The giant MPs 
present one or several broadband microbursts,which consist 
of narrowband (&v/v ~ 0.1) shots of a nanosecond duration. 
The giant IP's consist of the proportionally spaced narrow 
emission bands of microsecond lengths organized into several 
band sets. 

The radio profile components of pulsars outside of the 
MP used to be interpreted in terms of several geometrical 
models (e.g jHankins fc Cordeslll98ll ). which assumed a pe- 
culiar geometry of the pulsar (either orthogonality or align- 
ment of its magnetic and rotation axes) and placed the emis- 
sion regions of the components at different sites in the mag- 
netosphere. Recently this approa ch was complemented b y 
the concept of inward emission (|Dvks. Zhang fc Gill I2OO5I ) : 
although pulsar radiation is undoubtedly emitted along the 
magnetic lines, it may well be directed backwards, toward 
the neutron star. However the complicated radio emission 
pattern of the Crab pulsar remains beyond the framework 
of any plausible geometrical model. Moreover, the compre- 
hensive spectral, polarization and fluctuation data testify 
against the independent generation of the components and 
highlight the question of their physical nature. 

Recently we have proposed a physical mechanism of the 
PR and IP components in radio pulsars (|Petrovall2008al rbl). 
It is based on the propagation effects in the flow of pair 
plasma inside the open field line tube of a pulsar. The PR 
and IP are suggested to result from the induced scattering of 
the MP radiation into background. The scattered radiation 
strongly concentrates in the direction corresponding to the 
maximum scattering probability and forms separate compo- 
nents in the pulse profile. In the regime of a superstrong 
magnetic field, the scattered radiation is directed approxi- 
mately along the field, and, taking into account the effect of 
rotational aberration, one can recognize it as a PR compo- 
nent. In a moderately strong magnetic field, the radiation 
is predominantly scattered backwards - in the direction an- 
tiparallel to the velocity of the scattering particles, which 
stream along the magnetic lines, - and the scattered com- 
ponent can be identified with the IP. 

In application to the Crab pulsar, we have generalized 
the problem of induced scattering by taki ng into account the 
helical motion of the scattering particles jPetrova 2008c]). In 
the outer magnetosphere, the particles acquire gyration en- 
ergies as a result of resonant absorption of radio emission. 
They arrive at relativistic gyration at the very bottom of 



the resonance region, where the highest radio frequencies 
are absorbed. At higher altitudes, the character of the ra- 
dio wave scattering changes essentially: now the waves can 
be scattered between different harmonics of the particle gy- 
roferquency. It has been shown that the induced scattering 
between the states well below the resonance may still ac- 
count for the PR component in the radio profile of the Crab 
pulsar, whereas the scattering from the first harmonic to th e 
zeroth one may be responsible for the LFC (|Petrovall2008d ) . 

As the radiation of the PR and LFC appears well be- 
low the resonance, further on its way in the magnetosphere 
it is subject to the induced scattering in a moderately strong 
magnetic field. Similarly to the MP scattering, it may give 
rise to the additional backward components in the pulse pro- 
file. A preliminary analysis of the Crab polarization data 
suggests that the IP' results from the backward scattering 
of the PR, whereas the pair of HFCs is a consequence of the 
backscattering of the LFC. 

It should be noted, ho wever, that the simplified scat- 
tering geometry derived in iPetroval (|20 08b) cannot be di- 
rectly extended to the case of the Crab pulsar: neither of the 
backscattered components, including the IP, can be properly 
located in the pulse profile. Given that the backscattered 
radiation is directed along the magnetic lines and suffers ro- 
tational aberration, the scattered component should lag the 
original one by m ore than 180° in pulse longitude (see also 
iDvks et al.ll2005r ). At the same time, the IP of the Crab pul- 
sar is separated from the MP by only ~ 150°, the distance 
between the IP' and PR is almost the same and the separa- 
tion of the HFC 1,2 from the LFC, being larger than 180°, 
is still small to account for the scattering above the LFC 
origin. 

In the present paper, we overcome this seeming diffi- 
culty by means of a more accurate examination of the com- 
ponent formation. The paper presents a detailed study of the 
intensity transfer in the course of induced scattering of pul- 
sar radiation into the background. On this basis we elaborate 
our theory of the pulsar radio profile components outside of 
the MP, extending it to the case of the Crab pulsar, with 
an especial emphasis on the formation of its high-frequency 
components (the IP', HFC1 and HFC2). In Sect. Owe solve 
the problem of induced scattering of a narrow radio beam 
into a wide solid angle in the presence of a moderately strong 
magnetic field. With the angular distribution of the scat- 
tered intensity in hand, in Sect.|3]we simulate the profiles of 
the high-frequency components of the Crab pulsar, taking 
into account the effect of rotational aberration. The spec- 
tral, polarization and fluctuation properties of the scattered 
components are confronted with the observational evidence 
in Sect. [3] The results are briefly summarized in Sect. [5] 



2 INDUCED SCATTERING INTO A WIDE 
SOLID ANGLE 

2.1 Statement of the problem 

Pulsar magnetosphere contains the ultrarelativistic electron- 
positron plasma, which streams along the open magnetic 
lines. The radio emission originates deep in the magneto- 
sphere and subsequently propagates through the plasma 
flow. The radio wave scattering off the plasma particles 
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is one of the significant propagation effects in pulsars. As 
the brightness temperatures of pulsar radio emission are ex- 
tremely high, the indu ced scattering strongl y dominates the 
spontaneous one (e.g.. IWilson fc Reeslll978r ). 

The external magnetic field may affect the scattering 
process considerably provided that the frequency of inci- 
dent waves in the rest frame of the scattering particles is 
less than the particle gyrofrequency, u>' = 1^7 (1 — f3 cos 9) <C 
ujg = eB/mc (here u) is the wave frequency in the labo- 
ratory frame, (3 is the particle velocity in units of c, 7 is 
the Lorentz-factor of the particle motion, 7 = (1 — /3 2 )~ 1//2 , 
9 is the angle of incidence of the waves). The radio waves 
typically meet the cyclotron resonance condition, u)' = u>g, 
in the outer magnetosphere, and therefore the pulsar ra- 
dio emission is subject to both the magnetized and non- 
magnetized scattering. Th e non-magnetized indu c ed scat- 
tering was considered in IWilson fc Reel dl97Sl ); IWilson " 



in the laboratory frame the photon states involved in the 
scattering satisfy the condition 



017(1 — p cos 9) = 0117(1 — j3 cos 9i ). 



(1) 



ermg was considered m I W nson fc Kccs (iy/8J; Wilson 
19821 ; lHardv fc Melrose! (|l995h ; iLvubarsklT fc Petroval 



1996 ) ; iLvubarskv ( 2008[). and the magnetized s catter - 



ing was addressed in IBlandford fc Scharlemannl (| 19761) 



iLvubarskii fc Petroval (|l996l ); iPetroval l|2004al lbl. l2008al lbl). 

The presence of a strong ma gnetic field can modify the 
scatt ering process variously (e.g.. IBlandford fc Scharlemannl 
1 19761 ). In an extremely strong field, the perturbed motion of 
a particle in the fields of the incident wave is confined to 
the magnetic line. Due to such a character of the particle 
motion the scattering process has a number of peculiarities 
and is usually called the longitudinal scattering. In this case, 
only the waves of the ordinary polarization, with the electric 
vector in the plane of the ambient magnetic field, suffer the 
scattering, and the scattered waves also have the same polar- 
ization. Besides that, the waves are predominantly scattered 
forwards, at ~ I/7 to the magnetic field direction. In a mod- 
erately strong magnetic field, the components of the parti- 
cle's perturbed motion perpendicular to the magnetic line 
affect the characteristics of the scattered radiation signifi- 
cantly. Both the ordinary and extraordinary waves partici- 
pate in the scattering process, and the induced scattering is 
most probable in the backward direction. This is a so-called 
transverse scattering. The induced scattering switches be- 
tween the regimes on condition that 7 2 o/ /lu g — 1 (|Petroval 
2008b), both regimes being typical of the pulsar magneto- 
sphere. 

Pulsar radio emission is known to be highly directional. 
At any point of the pulsar emission cone it is concentrated 
into a narrow beam of the opening angle ~ I/7. Therefore 
the induce d scattering out of the beam dominates that inside 
the beam (jLvubarskii fc Petrova|[l996h . The former process 
may only hold if initially there is some background radiation 
outside of the beam. This background radiation may result 
from the spontaneous scattering out of the beam. Although 
the background intensity is extremely small as compared 
to the beam intensity, Ihg/Io = 10 -8 — 10 -12 , it may still 
trigger an efficient induced scattering of the beam radiation 
into the background. 

As soon as the induced scattering starts, the intensity 
of the background grows exponentially, ft is clear that a 
substantial growth may occur only for the states with the 
photon orientations #1 close to that corresponding to the 
maximum scattering probability, #™ ax . In the rest frame of 
the scattering particles, the frequencies of the incident and 
scattered photons are almost the same, u' ~ uj[. Hence, 



In our previous papers, we have considered the induced scat- 
tering between the two photon states satisfying equation ([1]), 
one of which represents the photon beam of a given fre- 
quency u> and orientation 9 (I/7 <C 9 < 1) and anot her one 
the sta te of the most probable scattering, 9i = Q™ 3 * (|Petroval 
2008a. be). Then the broadband picture of the scattering has 
been reconstructed from the two-state solutions at different 

id. 

It should be kept in mind, however, that the scattering 
efficiency depends not only on the orientation of the scat- 
tered photons, but also on the incident intensity, which is 
a function of frequency. For a given frequency of the back- 
ground photons, uji, different photon orientations 9\ corre- 
spond to different feeding frequencies u (cf. Eq. fT}) and, 
consequently, to different original intensities of the radio 
beam. If, say, #™ ax = n , smaller 9\ imply smaller u) and 
larger To, so that the scattering into these directions may 
be even more efficient. Thus, the location of the scattered 
components in the pulse profile can be affected by the radio 
beam spectrum. 

In the present paper, we take into account this effect 
and examine the induced scattering of a narrow radio beam 
of a fixed frequency into all conceivable states of the back- 
ground, which satisfy equation JTJ. As the beam photons of 
different frequencies are scattered independently, the exten- 
sion of our results to the broadband case is straightforward. 
The main motivation of our study is the necessity of mod- 
eling the profiles of the scattered components as applied to 
the Crab pulsar, which is known for its extremely steep ra- 
dio spectrum. We consider the transverse scattering, which 
gives rise to the backward components. As for the longitu- 
dinal scattering, its efficiency is a much stronger function of 
the orientation of the scattered photons, and the radio beam 
spectrum does not affect the profile location of the scattered 
components considerably. 



2.2 Intensity transfer 

Let us consider the transverse induced scattering in the 
frame corotating with the neutron star. The evolution of the 
photon occupation number n in the course of the s c atterin g 
is described by the kinetic equation (1) in IPetroval (l2008bl ). 
One can rewrite it in terms of intensity, i v — 2hv 3 n/c 2 , and 
integrate over the Lorentz-factor of the scattering particles, 
noticing that a detailed form of the particle distribution 
function does not play a crucial role and taking the mo- 
noenergetic distribution with some characteristic Lorentz- 
factor. The radio beam intensity is assumed to have the 
delta-functional angular distribution, and it is convenient to 
introduce the spectral intensity of the beam integrated over 
the solid angle, I v = J i„ bcam df2. Pulsar radio emission is 
generally a mixture of the two types of waves: the ordinary 
ones polarized in the plane of the ambient magnetic field 
(the A-polarization) and the extraordinary ones polarized 
perpendicularly to that plane (the B-polarization). For the 
sake of simplicity we assume that the radio beam consists of 
the waves of one polarization and interacts with both polar- 
ization states of the background. Then the intensity transfer 
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between the beam and the background can be presented as 
follows: 

^ = I »'(0i,M(<»s0-ra»0i)0 tf (0i,fr)dni I 



di J (6»i 



da; 



= » j (0i„ > 0i J^(cos 61 - cos 0i a )g tJ (8i a , 4>i J , 



dt ' ( ^'M = * J (6% , )/ l (cos 6>-cos ,^), (2) 



This infinite set contains the equations for the background 
photons of any conceivable orientations , 01,), / = 
a, f3,.... Here the intensities of the beam and the back- 
ground are normalized to the original radio beam intensity, 
I = I v /Iv and i = iv hs /Iu , da; = ii'^adr, 



4n e rg/3 v' 



(3) 



n e is the number density of the scattering particles, r e is 
the classical electron radius, 8 is the beam tilt to the mag- 
netic field, I/7 < 8 < 1, v' = 1/7(1 -/3cos0), i/ G = wg/2tt, 
r is the coordinate along the beam trajectory in the scat- 
tering region, the subscripts of the factor g 13 stand for the 
polarization states of the incident and scattered photons, 

/^(1 + sinW 1 -^ 2 ( 1 -^) 2 



/3 2 7V (3 2 j 4 vl 



sin 29 sin 29i cos A<f> 



ab . 2 A (1-777 ) 

<7 = (1 + cos A(p)- 



g = (1 + cos Acp)- 



P 2 ~t A nl 



g BB = 1 + sin 2 A(j>, 



(4) 



A<f> = 4>i — 4> is the azimuthal coordinate of the background 
photons with respect to that of the beam, r\ = 1 — j3 cos 6 and 
•qi = 1 — j3 cos 61 . In the set of equations ([2]) it is taken into 
account that initially the background intensity is very small 
as compared to the beam intensity and therefore the terms 
describing the induced scattering between the states of the 
background photons are omitted. One can see that the total 
intensity of the beam and the background is conserved, 



f + 



3 



i J dfil EE C » 1. 



(5) 



Furthermore, from the last two equations of the set ((2)1 one 
can find that for the two arbitrary orientations of the back- 
ground photons, (6i a , 4>i a ) an( i (#1,3 > 0is)j the intensities are 
related as 



,) 



i 3 o (6i ,c/>i 



.(6) 



Making the substitution 

i 3 (0i, , 01,) = ij (0i, , 0i, ) exp(y'g 13 (cos - cos $ h )) (7) 

and using equations ©-(O in any of the equations ([2]) for 
the background intensity, we obtain 



j=A,B 



Me 



y z g" 13 (cos 6— cos 6\ ) 



dfii 



(8) 



The solution y l (x) completely determines the intensity evo- 
lution in the course of the scattering. For any point of the 
scattering region, the beam intensity reads 



P = 



dx ' 



(9) 



whereas the angular distribution of the background intensity 
is given by equation (0. 

The differential equation |(HJ does not allow the analyt- 
ical treatment. The numerical solution is plotted in Fig. [TJi. 
Here and hereafter in this section we concentrate on the 
case of the A-polarized beam, keeping in mind that for any 
of the two polarizations the results are qualitatively similar 
and their quantitative difference is insignificant. The inten- 
sity evolution of the beam and the background is shown in 
Figs. [Tp and [TJ;, respectively. One can see that the inten- 
sity transfer becomes significant at the scattering efficiency 
x < 10. Note that for different directions 81 the background 
intensities arrive at the stage of saturation simultaneously, 
the final intensity values being substantially different. 

The difference of the intensity growth in the two polar- 
izations is introduced by the factor g lJ (see equation 0). 
As can be noticed from equation @, for 9i ~ tt the ratio 
g AA I g AB depends on 61 very weakly, but varies with Acj> 
within a factor of 4, which appears crucial for the relative 
intensity growth. Given that A(j> changes by 90° , the value of 
g AA I g AB is inverted. Thus, the dominant polarization of the 
scattered radiation is determined by Acj>. In Fig.[TJj, A<f) — 
is taken and the B-polarized intensities are shown, the A- 
polarized intensities being negligible. The final intensity of 
the background as a function of 81 is presented in Fig. [TJl. 
One can see that the angular distribution of the scattered 
intensity has a small but meaning width. 



3 APPLICATION TO THE CRAB PULSAR 

To get a notion about the efficiency of intensity transfer as 
a result of the transverse induced scattering in the Crab 
pulsar one can use the numerical estimate of the quantity 
3;(ee r) = /i°'ar, which is given by equation (12) in IPetroval 
(2008b). Substituting the parameters of the Crab pulsar (the 
period P — 0.033 s, the magnetic field strength at the surface 
of the neutron star _B* = 4 x 10 12 G, the radio luminosity at 
400 MHz L400 = 1.4 X 10 29 ergs -1 and the spectral index 
a = 3), taking the multiplic ity factor of the plasma k = 10 4 
|Harding fc Muslimovll20C)2r ) and the characteristic Lorentz- 
factor of the scattering particles 7 = 10 and assuming the 
angle of incidence of the radio emission 8 = 0.3, we find that 
at v — 1400 MHz x = 35. Although the characteristics of 
the pulsar plasma are very uncertain, one can conclude that 
in the Crab pulsar the scattering can indeed be efficient. 

The scattering efficiency depends on the frequency of 
the incident radiation. As can be seen from equation @, 
a x u 2 , and, taking into account the power-law spectrum 
of pulsar radiation, oc v~ a , we have x oc v 2 ~ a . Hence, 
if we are interested in the scattered component of a given 
frequency v\ (y\ = vr\jr\\), the scattering efficiency x oc 
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Vi~ a ■ Using this in equation Q, one can obtain the angular 
distribution of the scattered intensity at the frequency v\. 

It should be kept in mind, however, that the formalism 
of Sect.[2]is written in the frame corotating with the neutron 
star at the angular velocity ft. To analyse the observational 
consequences of the induced scattering it is necessary to per- 
form the transformation to the laboratory frame. As a result 
of this transformation, the corrections to all the quantities 
involved, except for the photon orientations, are of the sec- 
ond and higher orders in t/tl (here tl is the light cylinder 
radius) and can be ignored in our consideration. The effects 
of rotational aberration and retardation introduce the first- 
order corrections and both should be taken into account. 
Note that the aberration can cause not only quantitative 
but also qualitative changes in the profiles of the scattered 
components seen by an observer, and therefore this effect is 
of an especial interest. 



3.1 The effect of rotational aberration 

Generally, the rotational aberration acts to redistribute the 
intensity, shifting and concentrating the rays toward the 
direction of the rotational velocity f! x r. Note that the 
rays inclined at the angles 180° i( to ft x r are shifted 
in opposite directions, giving rise to the two humps on 
the intensity distribution. The effects of rotational aber- 
ration and retardation in pulsar magnetosph ere have been 

~1 

2001 • 



addressed in a number of papers (see, e.g.. iMorinJ 



Krishnamohan & Downs 1983; Hibschman & Arons 



Gangadhara fc Guptall200ll: iDvks. Rudak fc Hardindl20ol 



Wang et a 



1120061). Here w e proceed from the formalism de- 



Petroval (|2008bl ) for the rays backscattered along 



veloped in 

the magnetic lines, generalize that analysis for the rays of 
an arbitrary orientation and examine the influence of the 
rotational aberration on the intensity distribution. 

As the regime of transverse scattering holds in the outer 
magnetosphere, one can assume that the emission altitude 
is negligible as compared to the altitude of the scattering 
region and the ray geometry is dominated by the effects of 
magnetosphere rotation, which include the rotational aber- 
ration in the scattering region and the retardation because of 
ray propagation to the scattering point in the rotating mag- 
netic field. For the sake of simplicity we take that the mag- 
netic and rotational axes of the pulsar are perpendicular. 
The geometric scheme of the scattering is shown in Fig. [2] 
The ray fc emitted along the magnetic axis at t = at the 
point O comes to the point of scattering S at t = r/c, while 
the magnetic axis turns by the angle fir /c = r /tl - The polar 
angle of the point of scattering with respect to the instanta- 
neous magnetic axis is r/rz,, and, in the dipolar geometry, 
the ambient magnetic field vector b is inclined to the axis at 
the angle 3r/2rL- Then the angle between fc and b equals 
r/2rL- At t = r/c, the radiation is emitted at point O along 
the instantaneous magnetic axis and the scattered radiation 
arises at point S at the angle 3r/2r_L + 9i + A9 a to that di- 
rection (here A9 a is the correction for the rotational aberra- 
tion), so that the scattered radiation lags the original one by 
3r /2rL+9i+A9 a in rotational phase. It should be noted that 
the scattered component travels somewhat larger distance to 
the observer than the MP, Ar = r cos{n — 9\ — A8 a — r/2ri). 
Hence, the pulse longitude of the scattered component with 
respect to that of the MP reads 



A sc - Amp = ^— + 0i + A6 a - — cos(6»i + A9 a ). (10) 

2 TL Th 

Here we have retained only the quantities of the first order 
in r/r_L, since the higher-order effects are ignored in our 
consideration. 

The relativistic aberration is described by the well- 
known formula 



cos & 



COS £ c + f} r 



(11) 



1 + fi r COS £ c 

where and £ c are the wavevector tilts to the rotational ve- 
locity in the laboratory and corotating frames, respectively, 
and fi r is the rotational velocity in units of c, f3 r = |ft x r\/c. 
Taking into account that 



6 



2 



9l < I - 2F7: 



(12) 



and using equation (|lip . one can find the correction for the 
aberration A9 a = #i(Ci) — $i(£c). 

The results of numerical calculations based on equations 
(I10I) -(I12 [ ) are presented in Fig. [3] It shows the longitudinal 
location of the rays on the pulse profile versus their tilt 9i 
to the magnetic field b for different altitudes of the scat- 
tering region. One can see that the rays equidistant in 9± 
tend to concentrate in the two longitudinal regions of the 
profile, the effect being pronounced for large enough alti- 
tudes of the scattering region, in which case the rotational 
aberration is sufficiently strong. This trend may affect the 
resultant intensity distribution. The concentration of rays in 
the region closer to the MP may account for the emission 
bridge linking the MP and IP. Given that the second region 
of ray concentration does not coincide with the peak of the 
original intensity distribution, it may be responsible for a 
separate component in the pulse profile. 

The geometrical considerations leading to equation (|10|l 
concern the transverse scattering of the MP into the IP. 
However they can also be applicable to the backscattering of 
the PR and LFC, which are by themselves the results of the 
MP scattering in the longitudinal regime at lower altitudes 
in the magnetosphere. It can be shown that if the PR and 
LFC are subject to the transverse scattering far enough from 
their origin, r > ro, the scattering geometry is deterimned 
by the magnetosphere rotation rather than by the details of 
the previous scattering and the longitudinal separations of 
the resultant backscattered components from the PR/LFC 
locations on the pulse profile (Ai P / — Apr and Ahfc — Alpc) 
are still given by equation (|10p . 



3.2 Profiles of the backscattered components 

The numerically simulated profiles of the backscattered com- 
ponents of the Crab pulsar are shown in Fig. [4] The location 
and shape of the components are determined by the 3 pa- 
rameters: the scattering efficiency x c for a fixed 0i, the spec- 
tral index a of the incident radiation and the altitude of the 
scattering region, r/r_L. The scattering efficiency affects the 
peak location, since it determines the orientations of the rays 
arriving at the stage of saturation. Larger x c shift the scat- 
tered component toward later pulse longitudes. The steeper 
the spectral index of the incident radiation, the narrower is 
the scattered component. The quantity r /tl determines the 
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strength of the effect of rotational aberration, as is discussed 
above. 

As the LF C originates high in the magnetosphere 
(Petrova 2008c), its scattering to the HFC takes place close 
enough to the light cylinder, and the rotational aberration 
plays a significant role. As can be seen in Fig. it shifts 
the scattered component as a whole to later pulse longitudes 
(cf. Fig. [3} and gives rise to an additional peak (known as 
HFC2) in the intensity distribution. The PR is backscat- 
tered into the IP' at lower altitudes, where the rotational 
aberration is much less significant. The IP' peaks at ear- 
lier pulse longitudes (see Fig. [Jp) and, due to the extremely 
steep spectrum of the PR, «pr = 5, looks narrow. The scat- 
tering of the MP into the IP is qualitatively similar, but the 
MP needs much steeper spectrum than is observed to pro- 
vide the IP narrowness. It should be noted that the induced 
scattering may flatten the spectrum of the incident radia- 
tion, being more efficient at lower frequencies. Hence, the 
original spectrum of the MP may really be steeper. In our 
simulation presented in Fig. [4jj we have assumed omp = 4, 
but it still seems smaller that necessary. 

Note that for the scatterings PR^IP' and LFC^HFCs 
we assume the A-polarization of the incident radiation, in 
which case the scattered radiation is dominated by the B- 
polarization. The polarization type of the PR and LFC is un- 
ambiguously determined by the physics of longitudinal scat- 
tering which gives rise to these components (see Sect. [2~Tj) . 
The position angle shifts of the IP' and HFCs by 90° from 
those of their feeding component s is in line with the obser- 
vations (|Moffett fc Hankins|[l999l ). 



4 DISCUSSION 

We have suggested the scenario of the component formation 
in the Crab pulsar (see Fig. [S] for summary) . It is based 
on the induced scattering of the MP into the background, 
which takes place in the flow of the secondary plasma inside 
the open field line tube of the pulsar. Deep in the magneto- 
sphere, the magnetic field is so strong that the MP radiation 
is scattered in the longitudinal regime and gives rise to the 
component directed approximately along the streaming ve- 
locity of the scattering particles. It appears on the pulse 
profile ahead of the MP and can be identified with the PR 
(jPetro va 2008a. cj). As the magnetic field strength decreases 
with distance, at higher altitudes the scattering switches to 
the transverse regime and results in the component roughly 
antiparalle l to the particle velocity, which can be recognized 
as the IP (|Petrovall2008bl ). Resonant absorption of the ra- 
dio emission leads to the relativistic gyration of the scat- 
tering particles, and the scattering from the first harmonic 
of th e gyrofrequency to the zeroth one can account for the 
LFC (IPetrovall2008cr ). On the way in the magnetosphere, the 
PR and LFC are subject to the transverse scattering, which 
gives rise to the IP' and HFCs, respectively. 

According to equation JT]), the scattered radiation of 
a given frequency and different orientations is fed by the 
incident radiation of substantially different frequencies. The 
scattering efficiency x depends on the incident intensity, and 
hence, the steep spectrum implies an essential distinction 
of x at different pulse longitudes of the scattered compo- 
nent. It is the effect that proved to determine the peak lo- 



cations of the backscattered components of the Crab pul- 
sar (see Sect. [3j. In case of the MP scattering to the PR 
and LFC, however, the scattering efficiency appears much 
stronger function of the ray orientation, and the components 
always peak at 6\ ~ I/7, independently of the spectrum of 
the incident radiation. 

The observed distribution of the scattered intensity can 
be noticeably affected by the rotational aberration, which 
typically shifts and concentrates the rays toward later pulse 
longitudes. The effect is especially pronounced in the HFCs, 
which are believed to originate close to the light cylinder. In 
the framework of our model, they present a single compo- 
nent, with one of the peaks (the HFC1) being determined 
by the efficiency of induced scattering and another one by 
the rotational aberration. In the case of the MP scattering 
to the leading components (the PR and LFC), at a fixed 
altitude the scattered intensity has almost delta-functional 
angular distribution. Then the role of rotational aberration 
is to place these components ahead of the MP and to pro- 
vide their finite widths by shifting the radiation originating 
at somewhat different altitudes to different pulse longitudes. 

4.1 Frequency evolution of the profile 

The profile o f the Crab pulsar drastica lly changes with fre- 
quency (e.g., iMoffett fc Hankinslll996D : each of the compo- 
nents is visible only over the part of the radio frequency 
band. In our model, the frequency evolution of the compo- 
nents can be understood as follows. As the characteristic al- 
titudes of the scattering regions are related to the cyclotron 
resonance radius, lower frequencies are scattered somewhat 
higher in the magnetosphere. The scattering MP-^LFC 
holds at r ~ vl, and for the lowest frequencies the scatter- 
ing site is expected to lie beyond the light cylinder, in which 
case the stable component is not formed. Therefore the LFC 
appears on the pulse profile at high enough frequencies. The 
re-scattering of the LFC into the HFCs takes place inside 
the light cylinder only for still higher frequencies. Hence, at 
the highest radio frequencies the HFCs become prominent, 
whereas the LFC vanishes. 

The efficiency of the transverse scattering explicitly 
scales with frequency as x cx v 2 ~ a , and since qmp = 3, 
the scattering MP— >IP is believed to be significant only at 
low enough frequencies. At higher frequencies it is domi- 
nated by the scattering MP— >PR— >IP'. It should be pointed 
out that the latter process takes place at lower altitudes. 
Indeed, the longitudinal scattering (MP^PR) occurs lower 
than the transverse one (MP^IP), since the magnetic field 
strength decreases with distance. The frequency of the PR 
is much higher than that of the MP, i/pr ~ ^mp#mp/$pr ~ 
^mp#mp7 2 (see equation fTJ), and the PR tilt to the ambient 
magnetic field, Apr, rapidly increases along the trajectory 
from ~ I/7 to ~ r/tL, so that this component meets the 
condition of transverse scattering very soon. Thus, the IP' 
originates at lower altitudes than the IP, the former process 
being more efficient, in particular, because of larger number 
density of the scattering particles, n e oc r~ 3 . Probably, it is 
the altitudinal dependence of the scattering efficiency that 
provides stronger transverse scattering at higher frequencies 
(i.e. at lower altitudes) and accounts for the increasing spec- 
trum of the IP' and the flat spectrum of the HFCs. 

Note that the spectra of the components subject to the 
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scattering should be substantially altered by this process, 
so that the original spectra of the MP, PR and LFC differ 
from the observed ones and remain obscure. The spectral 
indices used in our simulations (see Fig. |4| are taken only 
for illustrative purposes. A detailed analysis of the spectra of 
the backscattered components is also beyond the framework 
of the present paper. 

The observations have revealed several fine effects in 
the high-frequency profiles: as the frequency increases, the 
IP' widens, the HFCs shift to later pulse longitudes, and 
the HFC2 becomes som ewhat stronger than the HFC1 
ijMoffett fc Hankinj[l999h . In principle, all these effects can 
be reproduced by means of a slight variation of the param- 
eters taken to construct the profiles in Fig. [4] However it 
should be kept in mind that our consideration is too simpli- 
fied, since we have assumed a fixed altitude of the formation 
of a backscattered component and neglected the width of an 
incident component as compared to its tilt to the ambient 
magnetic field. The fine effects should be considered in a 
more accurate model. 

4.2 Polarization properties and geometry of the 
pulsar 

The position angle swing across the profile of the Crab pul- 
sar cannot be explained in terms of the customary rotating 
vecto r model l|Moffett fc Hankindll999l ; iKarastergiou et al.l 
2004). Evidently, the assumption that all the profile com- 
ponents originate at the same altitude is not valid in this 
case. In the framework of our model of the Crab emission 
pattern, the profile components arise at different sites in the 
magnetosphere, the radiation can be directed antiparallel 
to the particle velocity, and the scattering can introduce a 
90°-shift of the position angle. Then the polarization picture 
roughly looks as follows. 

Very little change of the position angle across the MP 
testifies to the nearly central cut of the emission cone by the 
sight line. In the PR and IP components, which result from 
the MP scattering, the position angle is almost the same. 
The position angle of a scattered component reflects that of 
the incident radiation in the scattering region. Although the 
PR and IP originate at high enough altitudes, the position 
angle swing is still determined by the very small impact 
angle of the observer to the magnetic axis, /3, and is not 
modified substantially (i.e. (r/rz,) 2 sin a cos a <C f3, where a 
is the angle between the magnetic and rotation axes of the 
pulsar) . 

It should be noted that, in accordance with the physics 
of longitudinal scattering, the PR component should be 
characterized by complete linear polarization of the A-type. 
The ordinary polarization of the PR, along with the MP and 
IP, is confir med by the high-energy observa tions of the Crab 
PWN (e.g.. lLai. Chernoff fc Corded l200ll ). The transverse 
scattering of the PR leads to the IP' with almost complete 
B-polarization (cf. equation (f4)l). Apart from the 90°-shift, 
the position angle of the IP' is believed to roughly reflect 
that of the PR. 

The LFC shows an essentially distinct behaviour of 
the position angle: it increases linearly with pulse longi- 
tude, and the total shift from the MP position angle is 
Aip « Amp — Apr w 30° . This can be presumably attributed 
to the peculiar location of the LFC origin - close to the light 



cylinder. For an arbitrary angle a between the magnetic and 
rotational axes (given that (r/rz,) sin a cos a 3> 0) far from 
the emission region the rays make the angle (r/2ri,) sin a 
with the ambient magnetic field, and the position angle is 
Aip = (r/2rz,) cos a. Then the component scattered along 
the ambient magnetic field appears on the pulse profile at 
the longitude Alfc = Amp — (r/2rz,) sin a and has the same 
position angle, Aip = (r /2vl) cos a. (Let us note in passing 
that in an arbitrary geometry of the pulsar the backscat- 
tered components may also be visible provided that their 
orientation is determined by the pulsar spectrum rather 
than by the ambient magnetic field direction, as is sug- 
gested in the present paper.) Thus, the linear increase of 
the position angle with pulse longitude is indeed the case, 
Aip oc Amp — Alfc oc t/tl, but Aip ~ A A m 30° cannot be 
achie ved for any a. The value a « 60° cited in the literature 
fe.g.. lMoffett fc Hankinll 19991) yields only rough correspon- 
dence. In the more plausible case, a ~ 90°, we have Alp — > 0. 
The observed linear increase of Aip with A A oc r/r\t is sug- 
gestive of an other first-order effects, suc h as the action of the 
current flow (|Hibschman fc Arondl200ll ). though the higher- 
order effects close to the light cylinder cannot be excluded 
as well. 

The position angle of the HFCs is shifted by ~ 90° from 
that of the LFC, so that the A-polarized LFC is chiefly 
scattered into the B-polarization. In the two components, 
the HFC1 and HFC2, the position angle changes within the 
same range, confirming their origin as a single component 
affected by t he rotational aberrat i on. M ore observational de- 
tails given in lMoffett fc Hankind (| 1999T ) . e.g. different slopes 
of the position angle swing in the HFCs, cannot be com- 
pared to our model because of its restrictions (recall that 
we have assumed a fixed altitude of the component forma- 
tion and neglected the width of the component subject to 
scattering) . 

4.3 Giant pulses 

Similarly to the MP, all the backscattered components of 
the Crab pulsar - the IP, IP' and HFCs - exhi bit occasional 
giant pulses, whereas the PR and LFC do not (jCordes et al.l 
12004 Ijessner et alJlioolh . In the framework of our model, 
this can be understood as follows. As the efficiency of in- 
duced scattering is determined by the incident intensity, gi- 
ant MPs are scattered much stronger than the normal ones, 
and a substantial part of their energy may be transmitted 
to the scattered components, giving rise to the giant IPs, 
PRs and LFCs. Later on the giant PRs and LFCs are re- 
scattered into the IP's and HFCs, depositing most of their 
energy to these components and causing the giant IP's and 
HFCs. Thus, the scattering process transmits the energies of 
the giant MPs between different components, and these en- 
ergies are never stored in the intermediate components (the 
PR and LFC). It should be kept in mind that the energy of 
giant MPs is transferred between widely spaced frequencies 
(cf. equation (0), and, given that the giant MPs are suffi- 
ciently narrowband, at a given frequency the giant pulses in 
different profile components are not simultaneous. 

The induced scattering modifies the temporal and fre- 
quency structure of the incident radiation (see a lso iPetrova! 
2008b]). As is found in lHankins fc Eilekl (2007), the giant 
MPs and IP's of the Crab pulsar actually exhibit quite dis- 
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tinct structure. The giant MPs present the broadband mi- 
crobursts built of narrowband (5v/v ~ 0.1) nanoshots of a 
duration St ~ W/v ~ 10~ 8 - 10~ 9 s. The giant IP's con- 
sist of the proportionally spaced narrow emission bands of 
microsecond lengths, which drift noticeably toward higher 
frequencies. In our model, the distinction between the gi- 
ant MPs and IP's can be interpreted in terms of equation 
fl]) , which relates the frequency and orientation of the scat- 
tered radiation to those of the incident radiation. For the 
scattering MP— >PR— >IP' it takes the form 

v M p9 2 /2 = vpnh 2 = ^ip'(1 -/3cos6»i). (13) 

Differentiation of the extreme terms of equation (|13|l with re- 
spect to frequency yields A^mp / vmp = Afrp/ /vxpi . Thus, a 
nanoshot with 8v/u ~ 0.1 is scattered into the same normal- 
ized band 5v/u ~ 0.1, and the scattered radiation strongly 
concentrates toward the edge of the band because of the 
frequency dependence of the scattering efficiency. Then the 
groups of the nanoshots in the MP create the sets of emission 
bands in the IP'. 

The transformation of the angular scale of a nanoshot 
as a result of the scattering can be found by differen- 
tiating equation (I13|l with respect to angle: 2A9/9 — 
(3 sin 0i A0i / ( 1 - /3 cos Oi ) . Given that 7r - Ox -C 1, this yields 
A6»i » 4A6»/6>(tt - 0i) > AO. Thus, the angular scale of the 
scattered radiation increases by about two orders of mag- 
nitude. Keeping in mind that the radiation scattered from 
different nanoshots may overlap, one can expect that the 
group of nanoshots in the MP gives rise to the feature of 
microsecond length in the IP'. 

As can be seen from equation (|13|l . the radiation with 
somewhat different angles of incidence 9 (i.e. different parts 
of a nanoshot or the neighbouring nanoshots) is scattered 
to somewhat different frequencies f IP / , leading to the up- 
ward drift of the emission bands in the IP'. Thus, the pe- 
culiar temporal and frequency structure of the IP' can be 
attributed to the induced scattering of the MP. 



5 CONCLUSIONS 

We have examined the induced scattering of a bright nar- 
row radio beam into a wide solid angle. The radiation is 
scattered off the particles of the ultrarelativistic electron- 
positron plasma streaming along the open field lines in the 
magnetosphere of a pulsar. Our consideration is concerned 
with the transverse scattering in the presence of a moder- 
ately strong magnetic field. The intensity of the radio beam 
of a given frequency is transferred to the background radi- 
ation of the frequencies and orientations related by equa- 
tion fTJ, and the scattered intensity peaks in the direction 
antiparallel to the velocity of the scattering particles. The 
background radiation of a given frequency and different ori- 
entations 9i is fed by the beam radiation of different fre- 
quencies, and therefore the efficiency of intensity transfer to 
different 9i is determined not only by the scattering geom- 
etry but also by the corresponding original intensity of the 
radio beam. Hence, the angular distribution of the scattered 
radiation is affected by the radio beam spectrum, and for 
steep enough spectra of the pulsar radio emission it is the 
effect that determines the location of the scattered compo- 
nent on the pulse profile. 



Based on the numerical solution of the problem of inten- 
sity transfer we have simulated the profiles of the backscat- 
tered components of the Crab pulsar. It is suggested that 
the IP, IP' and HFCs result from the transverse scattering 
of the MP, PR and LFC, respectively. As the scattering re- 
gions lie in the outer magnetosphere, the observed profiles 
can be noticeably affected by the rotational aberration. A 
detailed treatment of this effect shows that not only the 
components are shifted in pulse longitude but also their in- 
tensities can be substantially redistributed. In particular, it 
is argued that the HFCs are a single component, with the 
second peak (the HFC2) arising as a result of the aberration. 

As the altitude of the scattering region is related to the 
radius of cyclotron resonance, higher frequencies are scat- 
tered somewhat lower in the magnetosphere. The number 
density of the scattering particles changes as r -3 , and the in- 
cident intensity as r -2 , so that at lower altitudes the scatter- 
ing is more efficient. The spectral behaviour of the backscat- 
tered components is presumably determined by the altitudi- 
nal dependence of the scattering efficiency: the IP' and HFCs 
appear in the high-frequency profiles and have flat or even 
increasing spectra. Note that at high frequencies the scatter- 
ing MP^IP is dominated by the scattering MP^PR-TP', 
which takes place at lower altitudes, whereas the scattering 
MP— >LFC^HFCs occurs inside the light cylinder and forms 
a stable component. 

Our model of the component formation in the Crab pul- 
sar is supported by the polarization data. The PR and the 
result of its backscattering, the IP', both have complete lin- 
ear polarization and the position angles differing by 90° . The 
position angles of the LFC and HFCs also differ by 90° . The 
flat position angle swing across the PR, MP and IP is sug- 
gestive of an almost central cut of the pulsar emission cone 
by the sight line. Given that the Crab pulsar is a nearly 
orthogonal rotator, the linear increase of the position angle 
with pulse longitude observed in the LFC as well as the to- 
tal position angle shift of ~ 30° from that of the MP can 
be attributed to the effect of the current flow close to the 
light cylinder. Note that all the scattered components can 
be observable for an arbitrary (but not very small) angle 
between the magnetic and rotation axes of the pulsar. 

In the framework of our model of the component for- 
mation in the Crab pulsar, the giant pulses characteristic 
of the MP can also be seen in the scattered components, 
since stronger incident intensities imply higher efficiency of 
intensity transfer to the scattered components. The lack of 
giant pulses in the PR and LFC can be understood as a 
consequence of their re-scattering to the IP' and HFCs. The 
induced scattering can also account for the observed differ- 
ence in the frequency and temporal structure of the giant 
MPs and IP's. 
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Figure 1. Intensity transfer as a result of induced scattering: a - numerical solution of the differential equation J5J; b - beam intensity 
in the A-polarization as a function of scattering efficiency; c - angular intensity of the background in the B-polarization vs. scattering 
efficiency for different directions 9\ (A<f> = 0); d - final angular intensity of the scattered radiation in the B-polarization vs. photon 
direction (A<f> = 0); 7 = 100, hg/Io = W~ 10 and 9 = 0.3. 
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Figure 2. Geometry of transverse scattering in pulsar magnctosphere. The pulsar is assumed to be an orthogonal rotator with the 
rotational axis perpendicular to the plane of the figure and the magnetic axis rotating counterclockwise. For more details see the text. 




Figure 3. Longitudinal location of the rays of equidistant orientations Q\ on the pulse profile with account for the rotational aberration 
at different altitudes in the magnctosphere 
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Figure 4. Simulated profiles of backscattered components: a - the HFC1 and HFC2 (the result of the scattering LFC A ^HFC S , where 
the superscripts mark the polarization type, the spectrum of the incident component is olfc = "3.15, the scattering efficiency at 
6»i = tt/2 is x c = 7, A<f> = 0, r/r L = 0.75); b - the IP' (the scattering PR A ^IP' B , ct PR = -5, x c = 10, A<f> = 0, r/r L = 0.45); c - the 
IP (the scattering MP B ^IP A , o M p = -4, x c = 5, A<f> = 0, r/r L = 0.6). 
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Figure 5. Summary of the component formation in the Crab pulsar. The marks L and T refer to the longitudinal and transverse 
of induced scattering, respectively. 



